Renewable  and  Sustainable  Energy  Reviews  28  (2013)  677-690 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Renewable  and  Sustainable  Energy  Reviews 

journal  homepage:  www.elsevier.com/locate/rser 


A  review  on  the  thermodynamic  optimisation  and  modelling 
of  the  solar  thermal  Brayton  cycle 

W.G.  Le  Rouxa,  T.  Bello-Ochende b'*,  J.P.  Meyer a 

a  Department  of  Mechanical  and  Aeronautical  Engineering,  University  of  Pretoria,  Private  Bag  X20  Hatfield,  Pretoria  0028,  South  Africa 
b  Department  of  Mechanical  Engineering,  University  of  Cape  Town,  Private  BagX3,  Rondebosch  7701,  South  Africa 


I  CrossMark 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  28  April  2013 
Received  in  revised  form 
31  July  2013 

Accepted  11  August  2013 
Available  online  6  September  2013 

Keywords: 

Solar 

Brayton 

Optimisation 

Entropy 

Modelling 


Many  studies  have  been  published  on  the  performance  and  optimisation  of  the  Brayton  cycle  and  solar 
thermal  Brayton  cycle  showing  the  potential,  merits  and  challenges  of  this  technology.  Solar  thermal  Brayton 
systems  have  potential  to  be  used  as  power  plants  in  many  sun-drenched  countries.  It  can  be  very 
competitive  in  terms  of  efficiency,  cost  and  environmental  impact.  When  designing  a  system  such  as  a 
recuperative  Brayton  cycle  there  is  always  a  compromise  between  allowing  effective  heat  transfer  and 
keeping  pressure  losses  in  components  small.  The  high  temperatures  required  in  especially  the  receiver  of 
the  system  present  a  challenge  in  terms  of  irreversibilities  due  to  heat  loss.  In  this  paper,  the  authors 
recommend  the  use  of  the  total  entropy  generation  minimisation  method.  This  method  can  be  applied  for 
the  modelling  of  a  system  and  can  serve  as  validation  when  compared  with  first-law  modelling.  The  authors 
review  various  modelling  perspectives  required  to  develop  an  objective  function  for  solar  thermal  power 
optimisation,  including  modelling  of  the  sun  as  an  exergy  source,  the  Gouy-Stodola  theorem  and  turbine 
modelling.  With  recommendations,  the  authors  of  this  paper  wish  to  clarify  and  simplify  the  optimisation 
and  modelling  of  the  solar  thermal  Brayton  cycle  for  future  work.  The  work  is  applicable  to  solar  thermal 
studies  in  general  but  focuses  on  the  small-scale  recuperated  solar  thermal  Brayton  cycle. 
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Nomenclature 

O 

Stefan-Boltzmann  constant,  W/m2  K4 

X 

period  time  for  one  oscillation,  s 

a 

constant,  - 

V 

volume,  m3 

A 

cross-sectional  area,  m2 

b 

constant,  - 

Subscripts 

BSR 

blade  speed  ratio,  - 

c 

constant,  - 

0 

surrounding/loss 

c 

specific  heat,  J/kg  K 

0 

zero  pressure  (ideal  gas)  for  cp 

COP 

coefficient  of  performance,  - 

1 

at  state  1 

cs 

gas  velocity,  m/s 

2 

at  state  2 

d 

degradation  constant,  - 

1-11 

refer  to  Fig.  1 

D 

turbine  diameter,  m 

a 

dead  (reference)  state 

e 

specific  exergy,  J/kg 

acc 

acceleration 

E 

exergy,  J 

b 

at  heat  transfer  boundary 

f 

sunlight  dilution  factor,  - 

B 

base 

Fd 

external  drag  force,  N 

BH 

Bahnke  and  Howard 

g 

gravity  constant,  m/s2 

be 

beam 

h 

specific  enthalpy,  J/kg 

c 

compressor 

I 

solar  irradiance,  W/m2 

col 

collector 

J 

polar  moment  of  inertia,  kg  m2 

cv 

control  volume 

k 

gas  constant  (cp/cv),  - 

d 

diffuse 

k 

thermal  conductivity,  W/m  K 

D 

destruction 

1 

thread  length,  m 

e 

exit 

L 

length,  m 

eff 

effective 

m 

mass  of  the  rotor  assembly,  kg 

eff 

efficiency 

m 

mass  flow  rate,  kg/s 

gen 

generation 

N 

speed,  l'pm 

i 

inlet 

N 

number  of  fins,  - 

j 

component  number 

NTU 

number  of  transfer  units,  - 

loss 

loss  into  the  environment 

P,P 

pressure,  Pa 

m 

mechanical 

Q 

heat  transfer  rate,  W 

max 

maximum 

r 

pressure  ratio,  - 

min 

minimum 

r 

distance  from  rotational  axis  to  thread's  attachment 

net 

net 

to  the  rotor,  m 

opt 

optimum 

R 

gas  constant,  J/kg  K 

P 

for  constant  pressure 

s 

specific  entropy,  J/kg  K 

P 

planet 

s 

solar  constant,  W/m2 

r 

receiver 

s 

entropy,  J/I< 

reg 

recuperator 

s 

entropy  rate,  W/I< 

rotor 

rotor 

t 

time,  s 

R 

reversible 

T 

temperature,  K 

s 

isentropic 

To 

environment  temperature,  I< 

s 

Sun 

U 

rotor  inlet  blade  tip  speed,  m/s 

sngl 

single 

V 

velocity,  m/s 

t 

turbine 

W 

power,  W 

V 

for  constant  volume 

X 

distance  in  x-direction,  m 

oo 

surrounding  area/free  stream 

y 

distance  in  y-direction,  m 

Z 

height,  m 

Superscripts 

Z 

constant,  - 

£ 

heat  exchanger  effectiveness,  - 

* 

>i 

efficiency,  - 

time  rate  of  change 

K 

dilution  factor  of  diffuse  radiation,  - 

CH 

X 

dimensionless  parameter  for  longitudinal  conduction, 

1.  Introduction 

Concentrated  solar  power  systems  use  the  concentrated  power 
of  the  sun,  as  a  heat  source  in  a  power  cycle,  to  generate  mechanical 
power.  Bejan  et  al.  [1]  implied  that  the  solar  heat  source  is  more 
suitable  than  the  isotope  and  nuclear  heat  sources  when  the  power 
plant  size  is  in  the  range  of  2-100  kW.  Solar  thermal  power  cycles 
have  potential  to  be  used  in  many  sun-drenched  countries.  One  of 


these  cycles  is  the  solar  thermal  Brayton  cycle.  This  cycle  can  be 
very  competitive  in  terms  of  efficiency,  cost  and  environmental 
impact.  Chen  et  al.  [2]  showed  that  the  Brayton  cycle  is  definitely 
worth  studying  when  comparing  its  efficiency  with  those  of  other 
power  cycles.  Mills  [3]  predicted  that  emphasis  may  shift  from 
Dish-Stirling  technology  to  solarised  Brayton  micro-turbines  due  to 
lower  Brayton  costs,  as  a  result  of  high  production  quantities  in  the 
current  market.  An  open  solar  thermal  Brayton  cycle  uses  air  as 
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working  fluid,  which  makes  this  cycle  very  attractive  for  use  in 
water-scarce  countries. 

In  a  solar  thermal  Brayton  cycle,  a  black  receiver  is  usually 
mounted  at  the  focus  of  a  sun-tracking  parabolic  dish  reflector  or 
mirror  field,  so  that  the  reflected  beam  irradiance  can  be  absorbed 
and  converted  into  heat  [4],  A  gas  is  heated  directly  or  indirectly 
from  the  concentrated  solar  radiation  and  it  is  then  used  in  the 
power  cycle  to  power  a  turbine.  A  major  disadvantage  of  the  solar 
thermal  Brayton  cycle  is  the  high  receiver  operating  temperatures 
required  to  get  reasonable  efficiencies.  Most  Brayton  cycles  are 
not  self-sustaining  at  operating  temperatures  below  480  °C  [4], 
Narendra  et  al.  [5]  noted  that  to  make  solar  thermal  power  a  real 
success  it  is  important  that  extensive  work  be  carried  out  in  the 
field  of  receiver  materials  and  their  constraints.  The  material 
constraints  of  the  solar  receiver  contribute  extensively  to  exergy 
loss. 

A  recuperated  solar  thermal  Brayton  cycle  allows  for  lower 
compressor  pressure  ratios  and  higher  efficiency.  A  recuperator 
(heat  exchanger)  is  used  to  extract  heat  from  the  turbine  outlet 
and  transfer  it  to  the  cold  stream  before  it  is  heated  in  the  solar 
receiver.  The  highest-efficiency  Brayton  cycles  are  recuperative 
cycles  with  low  compressor  pressure  ratios.  If  recuperation  is  not 
used,  high  compressor  pressure  ratios  are  required  to  provide  high 
efficiency  [4], 

Different  configurations  of  the  solar  thermal  Brayton  cycle  can 
be  found  in  the  literature.  The  direct  solar  thermal  Brayton  cycle, 
operating  at  relatively  low  pressure,  requires  a  large,  hot  gas  solar 
receiver  due  to  low  gas  heat  transfer  coefficients  [6],  An  indirect 
system  uses  an  extra  heat  exchanger  to  extract  heat  from  a  fluid 
system  in  direct  contact  with  the  solar  radiation  at  the  solar 
receiver.  This  configuration  allows  for  smaller  solar  receivers  to  be 
used.  In  an  open  Brayton  cycle,  air  as  working  fluid  is  free  to  enter 
and  leave  the  system.  For  open  Brayton  cycles,  it  is  crucial  to  have 
natural  air  movement  past  the  system  to  prohibit  the  reinjection 
of  the  warm  exhaust  gas  into  the  system  [4],  A  radiator  is  used  in  a 
closed  Brayton  cycle.  Helium  has  been  proposed  for  closed  solar 
Brayton  cycles  because  of  its  high  heat  transfer  capability  and 
because  it  is  inert. 

Heat  exchangers  can  also  be  used  for  intercooling  and  reheat¬ 
ing  in  a  solar  thermal  Brayton  cycle.  With  multi-staging  a  number 
of  compression  and  expansion  stages  are  combined  in  series  with 
coolers  and  heaters  respectively  [4],  According  to  Stine  and 
Harrigan  [4],  the  major  advantage  of  multi-staging  is  that  the 
cycle  can  have  the  high  efficiency  associated  with  a  recuperated 
low-pressure  ratio  cycle,  without  having  the  extremely  large 
recuperator  as  required  for  a  single-stage  cycle  of  the  same  power 
output.  Maximum  efficiency  is  achieved  when  equal  pressure 
ratios  are  maintained  across  each  compressor  and  each  turbine 
stage.  The  Ericsson  cycle  has  the  potential  of  attaining  Carnot 
efficiency  when  a  recuperator  is  used  [4],  Thus,  the  simplest  form  a 
solar  thermal  Brayton  cycle  is  an  open  system  consisting  of  a 
compressor,  turbine  and  solar  receiver.  More  complex  forms  of  this 
cycle  can  include  one  or  more  recuperators,  intercoolers,  rehea¬ 
ters,  radiators,  indirect  receiver  cycles  or  all  these. 

According  to  Stine  and  Harrigan  [4],  there  are  three  major 
losses  in  actual  Brayton  cycle  engines:  duct  pressure  losses,  turbo¬ 
machine  efficiencies  and  recuperator  effectiveness.  In  the  solar 
thermal  Brayton  cycle  major  heat  losses  also  occur  at  the  receiver, 
recuperator  and  air  exhaust  or  radiator.  No  matter  what  specific 
configuration  of  the  solar  thermal  Brayton  cycle,  the  influencing 
factors  to  the  performance  thereof  translate  mostly  into  the 
irreversibilities  of  the  system,  which  can  be  minimised.  Many 
studies  have  been  published  on  individual  components  used  in  a 
typical  solar  thermal  Brayton  cycle  and  their  optimisation, 
although,  often  these  components  are  optimised  individually 
and  not  as  a  part  of  a  whole  system  with  system  performance  as 


objective  function.  In  this  paper,  new  perspectives  on  the  optimi¬ 
sation  and  modelling  of  the  solar  thermal  Brayton  cycle  in  terms 
of  the  second  law  of  thermodynamics  are  given,  conflicts  in  the 
literature  are  shown  and  recommendations  for  future  work  are 
given. 

2.  Optimisation  of  the  solar  thermal  Brayton  cycle 

2.3.  The  second  law  of  thermodynamics 

Eq.  (1)  shows  an  exergy  balance  for  a  control  volume  [1,7], 
The  definitions  for  specific  exergy  transfers  at  inlets  and  outlets 
are  given  in  Eq.  (2)  [1] 

^=l(^-^Qj-(wcv-p0^j  +  'Z’niei-'Zmeee-ED  (1) 

e  =  (h-h0)-T0(s-s0)+1/2V2+gz+eCH  (2) 

There  is  a  fundamental  difference  between  energy  and  exergy: 
exergy,  unlike  energy,  uses  the  environment  as  its  reference  point. 
The  exergetic  approach  is  useful,  since  one  would  like  to  know 
what  the  optimal  possibilities  are  for  us  as  inhabitants  of  this 
environment  instead  of  inhabitants  of  infinite  space. 

Entropy  generation  is  a  path  dependant  property  and  should 
not  be  confused  with  the  thermodynamic  property  entropy  change 
(Sj-S,)  [1],  Eq.  (3)  [7]  gives  the  expression  for  the  balance  of 
entropy  for  a  control  volume 

miSi-ZmeSe+Z^+Sgen  (3) 

Consider  the  design  of  lifting  a  weight  with  a  pulley  system: 
entropy  generation  is  proportional  to  additional  work  wasted  for 
a  bad  pulley  system  design  [8],  Entropy  generation  can  be  mini¬ 
mised.  For  the  entropy  change  of  an  ideal  gas,  Eq.  (4)  [7]  can  be 
used  with  constant  specific  heat 

T  P 

s2  -S]  =  cp0  In  ~R  In  ^  (4) 

1 1  n 

2.2.  Exergy  analysis 

Exergy  analysis  is  often  applied  to  measure  the  performance  of 
a  solar  system.  Izquierdo  et  al.  [9]  for  example,  studied  the  entropy 
generated,  the  exergy  destroyed  and  the  exergetic  efficiency  of 
lithium-bromide  absorption  thermal  compressors  driven  by  solar 
thermal  collector  heat.  Total  entropy  generation  rate  for  the 
system  was  found  by  summing  the  entropy  generation  rates  from 
various  entropy  generation  mechanisms  in  the  system.  Narendra 
et  al.  [5]  presented  an  exergetic  analysis  of  a  solar  thermal  power 
system  using  the  Rankine  cycle.  This  analysis  was  done  for  the 
whole  system  to  show  the  irreversibilities  at  each  part  of  the 
system.  The  collector-receiver  assembly  was  found  to  be  the  part 
with  the  highest  losses.  Barrett  and  Reid  [10]  showed  that  the 
compressor  and  recuperator  are  the  main  contributors  to  the  total 
entropy  generation  rate  in  a  100  kWe  Closed-Brayton-Cycle  Space 
Power  System  optimised  for  minimum  mass.  Jubeh  [11]  did  an 
exergy  analysis  for  an  open  regenerative  Brayton  cycle  with 
isothermal  heat  addition  and  an  isentropic  compressor  and  tur¬ 
bine.  The  effect  of  two  heat  additions  in  a  gas  turbine  engine, 
rather  than  one,  was  analysed  from  the  view  point  of  the  second 
law  of  thermodynamics.  It  was  found  that  the  effect  of  isothermal 
heat  addition  on  the  performance  of  a  recuperated  Brayton  cycle  is 
positive  for  low  compressor  pressure  ratios,  low  ambient  tem¬ 
peratures,  and  high  turbine  inlet  temperatures.  Jubeh  [11]  empha¬ 
sised  that  first-  and  second-law  analyses  together  are  crucial  for 
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understanding  and  explaining  the  effect  of  any  parameter  on  the 
performance  of  a  thermal  system. 

2.3.  Entropy  generation  minimisation 

Thermodynamic  optimisation  or  entropy  generation  minimisa¬ 
tion  (EGM)  can  be  applied  in  an  exergy  analysis  to  optimise  the 
performance  of  a  system.  To  maximise  the  power  output  of  a  power 
plant  is  equivalent  to  minimising  the  total  entropy  generation  rate 
associated  with  the  power  plant  [12,13,14-16].  According  to  Bejan 
[12],  to  minimise  the  irreversibility  of  a  proposed  configuration,  the 
analyst  must  understand  the  relationship  between  temperature 
differences  and  heat  transfer  rates,  and  between  pressure  differ¬ 
ences  and  mass  flow  rates.  The  entropy  generation  rate  or  thermo¬ 
dynamic  non-ideality  of  the  design,  Sgen,  must  be  described  in  terms 
of  the  topology  and  physical  characteristics  of  the  system,  such  as 
dimensions,  speeds,  shapes  and  materials  [12], 

According  to  Bejan  [17],  EGM  experienced  astounding  growth 
during  the  1980s  and  1990s  in  both  engineering  and  physics  fields. 
The  EGM  method  relies  on  the  simultaneous  application  of  heat 
transfer  and  engineering  thermodynamics.  It  account  for  the 
inherent  irreversibility  due  to  heat,  mass  and  fluid  flow  processes 
of  devices  and  installations  [17], 

EGM  has  been  applied  in  various  fields.  The  first  power 
generation  field  to  use  EGM  simulations  regularly  was  that  of 
solar  power  plants.  It  was  found  that  an  optimum  coupling 
between  a  solar  receiver  and  power  cycle  exists,  so  that  the  power 
output  is  a  maximum  [17].  Shiba  and  Bejan  [18]  optimised  a 
counterflow  heat  exchanger  that  served  as  a  condenser  in  a 
vapour-compression-cycle  refrigeration  system  for  environmental 
control  of  aircrafts.  It  was  shown  that  the  minimisation  of  the  total 
power  requirement  was  completely  equivalent  to  the  minimisa¬ 
tion  of  the  entropy  generation  rate  in  the  entire  system.  According 
to  Torres-Reyes  et  al.  [19]  thermodynamic  optimisation  (or  EGM) 
can  be  applied  to  the  thermal  design  of  solar  heating  systems.  In 
their  paper,  the  architecture  of  a  solar  collector  with  a  specific 
geometry  and  specific  arrangement  between  the  absorber  surface 
and  the  fluid  duct  was  determined  based  on  global  optimisation. 
Using  the  theory  of  entropy  generation  minimisation,  Chen  et  al. 
[20]  optimised  the  internal  and  external  irreversible  closed  recup¬ 
erative  Brayton  heat  pump  cycle  coupled  to  constant  and  variable 
temperature  heat  reservoirs.  They  found  that  there  exists  an 
optimum  matching  between  the  effectivenesses  of  the  hot  side 
and  cold  side  heat  exchangers  of  the  system  and  the  regenerator 
and  that  for  maximum  COP  of  the  heat  pump,  an  optimum 
pressure  ratio  exists.  Zheng  et  al.  [21]  established  an  irreversible 
cycle  model  of  a  solar-driven  Braysson  heat  engine  (a  hybrid 
power  cycle  based  on  the  conventional  Brayton  cycle).  They  found 
that  for  maximum  efficiency,  there  is  an  optimum  value  of  the 
solar  collector  operating  temperature  for  a  given  set  of  operating 
parameters.  Wu  et  al.  [22]  also  presented  such  a  system  in  which 
the  temperature  dependent  heat  capacity  of  the  working  fluid,  the 
radiation-convection  heat  losses  from  the  solar  collector  and  the 
irreversibilities  due  to  heat  transfer  and  non-isentropic  compres¬ 
sion  and  expansion  processes  were  taken  into  account.  They  found 
that  the  overall  maximum  efficiency  and  optimum  operating 
temperature  of  the  collector  decrease  as  the  heat  losses  of  the 
collector  increases. 

According  to  Salamon  et  al.  [23],  minimising  the  rate  of  entropy 
generation  and  maximising  the  power  both  push  an  operation 
toward  minimum  wastefulness,  while  minimum  entropy  genera¬ 
tion  and  maximum  power  output  are  opposite  when  considering 
frugality.  Therefore,  minimum  entropy  generation  is  considered  as 
the  objective  of  the  conservationist  and  maximum  power  as  the 
objective  of  the  industrialist  when  for  example,  managing  an 
existing  power  plant  or  driving  a  motor  vehicle  [23],  In  this  paper, 


however,  the  term  entropy  generation  minimisation  refers  rather, 
to  the  optimisation  of  geometries  of  a  specific  design  such  that  a 
system  performs  optimally.  For  example,  the  geometries  of  the 
power  producing  parts  of  the  engine  of  a  motor  vehicle  can  be 
designed  in  such  a  way  that  it  produces  maximum  power  output 
when  a  specific  amount  of  fuel  is  available. 


2.4.  Optimisation  results  and  influencing  factors 

The  Brayton  cycle  and  its  optimisation  have  often  been  inves¬ 
tigated  in  the  literature.  Chen  et  al.  [24]  gave  the  power  output 
and  the  efficiency  of  a  regenerative  and  closed  Brayton  cycle  as 
functions  of  pressure  ratios,  reservoir  temperatures,  heat  exchan¬ 
ger  effectiveness,  compressor  and  turbine  efficiencies,  and  work¬ 
ing  fluid  thermal  capacitance  rates.  They  found  that  the  maximum 
power  output  of  the  system  is  strongly  dependant  on  the  effec¬ 
tiveness  of  the  recuperator  and  that  the  maximum  power  output  is 
attained  when  the  system  hot  and  cold  side  effectivenesses  are 
arranged  in  a  particular  way.  Wu  et  al.  [25]  studied  the  perfor¬ 
mance  of  an  endoreversible  regenerative  Brayton  heat  engine  with 
focus  on  minimising  irreversibilities  at  the  hot  and  cold  side  heat 
exchangers  of  the  system  (for  example,  the  solar  receiver  and  the 
radiator)  and  the  recuperator.  They  found  maximum  net  power 
output  of  the  system  when  the  hot  and  cold  side  heat  exchangers 
and  the  recuperator  are  ideal.  However,  these  conditions  require 
infinitely  large  components,  making  the  system  impractical. 

Stevens  et  al.  [26]  and  Stevens  and  Baelmans  [27]  stated  that  a 
recuperator  is  often  used  to  improve  the  overall  cycle  efficiency  of 
a  standard  gas  turbine  and  that  a  high  heat  exchanger  effective¬ 
ness  and  low  pressure  drop  are  required  to  achieve  maximum 
cycle  efficiency.  According  to  these  authors,  finding  a  compromise 
between  these  conflicting  requirements  is  the  main  challenge  in 
recuperator  design.  Stevens  and  Baelmans  [27]  found  that  for 
maximum  cycle  efficiency,  the  cold  and  hot  side  pressure  drops  of 
a  micro-recuperator  used  in  a  gas  turbine  cycle  are  uniquely 
connected,  since  their  ratio  depends  primarily  on  the  compressor 
pressure  ratio.  Stevens  et  al.  [26]  also  found  that  the  recuperator 
effectiveness  should  be  as  high  as  possible  and  that  the  pressure 
drop  should  be  preferably  located  at  the  cold  side  of  the  recup¬ 
erator.  Thus  they  concluded  that  the  hot  channels  should  be  larger 
than  the  cold  channels.  They  also  concluded  that  for  a  fixed 
recuperator  volume  optimisation,  the  recuperator  should  be  as 
short  as  possible  with  as  high  as  possible  cross-sectional  area. 

By  comparing  and  considering  three  different  gases,  Riazi  and 
Ahmed  [28]  presented  an  analysis  on  the  effect  of  the  ratio  of 
specific  heats  ( k )  of  a  closed  regenerative  solar  Brayton  cycle.  They 
found  that  a  fluid  with  a  high  ratio  of  specific  heats,  for  example 
air  (k=1.4),  has  peak  efficiencies  at  lower  compressor  pressure 
ratios  while  a  fluid  with  a  low  ratio  of  specific  heats,  for  example 
CF4  (fc=1.14),  has  peak  efficiencies  at  higher  compressor  pressure 
ratios.  Also,  a  fluid  with  a  high  ratio  of  specific  heats  performs  well 
when  the  temperature  ratio  (ratio  between  highest  and  lowest 
temperature  in  cycle)  is  also  high,  while  it  performs  poor  when  the 
temperature  ratio  is  low.  In  comparison  with  air,  a  fluid  with  a  low 
ratio  of  specific  heats  performs  better  at  low  temperature  ratios. 

Roco  et  al.  [29]  presented  a  general  theoretical  framework  for  a 
recuperative  gas  turbine  cycle  with  external  and  internal  irrever¬ 
sibilities,  in  order  to  study  the  maximum  power  output,  maximum 
efficiency,  efficiency  at  maximum  power,  power  at  maximum 
efficiency,  and  optimum  pressure  ratios.  The  optimal  operating 
conditions  were  investigated  in  terms  of  the  isentropic  efficiencies 
of  the  compressor  and  turbine,  pressure  drops  in  the  hot  and  cold 
streams,  and  effectiveness  of  the  heat  exchangers.  They  found  that 
the  maximum  efficiency  and  maximum  power  operating  points 
are  coincident  at  a  recuperator  efficiency  of  50%. 


W.G.  Le  Roux  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  28  (2013)  677-690 


681 


Cheng  and  Chen  [30]  used  thermodynamic  optimisation  to  opti¬ 
mise  the  power  output  of  an  endoreversible  intercooled  Brayton  cycle 
coupled  to  two  heat  reservoirs  with  infinite  thermal  capacitance  rates. 
The  effects  of  intercooling  on  the  maximum  power  of  an  endorever¬ 
sible  Brayton  cycle  were  examined.  They  found  that  an  endoreversible 
intercooled  Brayton  cycle  is  better  than  an  endoreversible  simple 
Brayton  cycle  without  lowering  the  thermal  efficiency.  A  similar  work 
was  also  done  by  Wang  et  al.  [31  ].  They  found  that  there  exists  an 
optimal  intercooling  pressure  ratio  for  maximum  power. 

Zhang  et  al.  [32]  established  a  model  in  which  the  heat  losses  of 
the  solar  collector  and  the  external  and  internal  irreversibilities  of 
a  solar-driven  Brayton  heat  engine  were  taken  into  account.  Thus, 
this  model  was  done  for  a  solar  thermal  Brayton  cycle  without  a 
recuperator.  It  was  found  that  the  performance  characteristics  of 
an  irreversible  solar  thermal  Brayton  heat  engine  are  different 
from  those  of  a  solar-driven  heat  engine  consisting  of  an  endor¬ 
eversible  Brayton  heat  engine.  It  was  found  that  the  larger  the  heat 
losses  in  the  solar  collector,  the  lower  the  optimum  operating 
temperature  of  the  solar  collector  and  the  smaller  the  maximum 
overall  efficiency  of  the  system. 


2.5.  Optimisation  of  the  solar  thermal  Brayton  cycle  using  EGM  and 
geometry  optimisation 

In  most  of  the  abovementioned  analyses  using  EGM  (or  mini¬ 
misation  of  irreversibilities),  however,  geometry  optimisation  was 
not  done  (except  for  example  [26]).  The  effects  of  component 
geometries  were  thus  not  included.  The  many  influencing  factors 
to  the  performance  of  the  system  can  be  much  better  understood 
with  the  use  of  entropy  generation  minimisation  and  the  optimi¬ 
sation  of  geometries  for  each  unique  case  of  the  solar  thermal 
Brayton  cycle. 

Geometric  optimisations  for  the  solar  thermal  Brayton  cycle 
were  done  by  Le  Roux  et  al.  [33-35],  Le  Roux  et  al.  [33]  found  that 
the  optimised  systems’  irreversibilities  were  spread  throughout 
the  system  in  such  a  way  that  the  internal  irreversibility  rate  was 
almost  three  times  the  external  irreversibility  rate.  The  geometries 
of  a  modified  cavity  receiver  and  plate-type  recuperator  in  an  open 
and  direct  solar  thermal  Brayton  cycle  were  optimised  using  EGM 
and  is  shown  in  the  following  example. 

When  taking  a  control  volume  around  an  open  and  direct  solar 
thermal  Brayton  system  and  assuming  a  steady  state,  it  can  be 
determined  where  exergy  is  crossing  the  boundary.  This  is  shown 
in  Fig.  1  [33-35].  An  exergy  analysis  was  performed  for  the  system. 
Take  note  that  the  inlet  and  outlet  temperatures  of  the  system  are 
not  the  same,  but  the  pressures  are.  Also  note  that  Tw=Tn  as  it  is 
assumed  that  the  control  volume  boundary  is  very  close  to  the  hot 


Fig.  1.  Control  volume  around  the  open  and  direct  solar  thermal  Brayton  cycle  to 
perform  an  exergy  analysis  [33-35], 


stream  exit  of  the  recuperator.  The  following  equations  arise 
which  can  be  used  as  objective  functions  [33-35]: 

T0Sgen  =  mermee+  (l  Q-  (l  -Q  QlossJ- Wt  +  Wc 

:.  T0Sgen  =  m  \Lhi-he-T0(Si-Se)  +  ](Vf-V^)+g(Zi-Ze)j 

+  (1-p)Q*-(l-QatoMJ-Wnet  (5) 

or 

Wnet  =  Wt-Wc  =  -T0Sgen+^(Vf-V^)  +  mg(Zi-ze) 

+  ^-^^-(l-^Qtossj  +  mCpodj-Tnl-riiToCpo  ln(-P-) 

(6) 

Furthermore,  the  total  internal  entropy  generation  rate  of  the 
system  can  be  written  in  terms  of  the  sum  of  the  entropy  generation 
rate  of  each  component  and  duct  in  the  system.  The  entropy 
generation  rate  of  each  component  is  in  turn  a  function  of  geometry 
variables.  Thus,  we  have  the  total  internal  entropy  generation  rate 
as  [33-35] 

Sgen  =  —  Sgen(compressor)  +  $ gerUductZi)  +  Sgen(recuperator)  +  Sgen(duct45) 

T  Sgen(receiver )  3“  Sgen(duct67 )  T  $ gen{turbine )  T  $gen(duct 89)  T  SgeMQthetj 

(7) 

Eq.  (7)  can  then  be  substituted  into  Eq.  (6)  to  get  the  net  power 
output  for  the  open  and  direct  solar  thermal  Brayton  cycle.  The  net 
power  output  is  then  written  in  terms  of  the  total  entropy  genera¬ 
tion  rate  of  each  of  the  components  and  ducts  in  the  system.  Note 
the  cancelation  of  the  terms  (T0/Tbj)Qi0ssj  and  (To/T*)Q*  in  the 
objective  function  when  the  entropy  generation  terms  from  Section 
3  are  added  into  Eqs.  (5)— (7).  This  equation  can  then  be  written  in 
terms  of  all  the  geometiy  variables  in  the  system.  This  equation  for 
the  net  power  output  is  the  objective  function,  which  should  be 
maximised  by  optimising  the  geometiy  variables  that  describe 
the  temperatures  and  pressures  at  each  point  in  the  system,  subject 
to  global  constraints.  The  components  of  Eq.  (7)  are  discussed  in 
Section  3. 

Eq.  (6)  is  also  found  in  a  similar  form  in  various  exergy  analyses 
of  existing  systems  for  example  by  Karsli  [36].  A  similar  approach 
for  a  nuclear  power  plant  was  taken  by  Gomez  et  al.  [37]  so  that 
the  total  entropy  generation  rate  was  written  in  terms  of  the 
power  output  of  the  system.  A  similar  exergy  analysis  was  also 
done  by  Jubeh  [  11  ].  In  these  systems,  however,  exergy  analysis  was 
done  without  optimisation. 

In  the  following  section  the  building  blocks  for  structuring  a 
unique  objective  function  for  the  optimisation  of  the  component 
geometries  of  a  specific  solar  thermal  Brayton  cycle  is  given.  The 
many  influencing  factors  to  the  performance  of  a  specific  solar 
thermal  Brayton  cycle  with  unique  components  can  be  much  better 
understood  with  the  use  of  entropy  generation  minimisation. 

3.  Structuring  the  objective  function  for  solar  thermal  Brayton 
cycle  optimisation 

In  this  section,  the  importance  of  modelling  different  sections  in 
a  system  so  that  an  objective  function  can  be  derived  is  emphasised. 
The  modelling  of  these  various  sections  and  components  are 
discussed  and  conflicts  and  new  perspectives  are  highlighted. 

3.1.  The  Gouy-Stodola  theorem 

The  Gouy-Stodola  theorem  (Eq.  (8))  states  that  the  lost  avail¬ 
able  work  (see  Eq.  (1))  is  directly  proportional  to  the  entropy 


682 


W.G.  Le  Roux  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  28  (2013)  677-690 


generation  in  a  system  [1,8,38],  This  theorem  allows  an  analyst  to 
describe  a  quantity  to  be  maximised,  such  as  the  power  output  of  a 
system,  in  terms  of  the  total  entropy  generation  rate  in  a  system. 
The  Gouy-Stodola  theorem  is  named  after  Gouy  (1889)  and 
Stodola  (1910)  who  recognised  its  importance.  The  Gouy-Stodola 
theorem  has  been  applied  and  referred  to  in  various  solar  thermal 
studies  including  for  example  Cervantes  and  Torres-Reyes  [39], 
Amhalhel  and  Furmariski  [40],  Sieniutycz  and  Kuran  [41]  and 
Hepbasli  [42] 

Ed  —  To  Sgen  (8) 

According  to  Lampinen  and  Wiksten  [43],  in  all  previous  studies, 
as  well  as  in  the  theories  given  by  Bejan  [44]  and  by  Tolman  and 
Fine  [45],  a  general  method  of  finding  the  correct  temperature  to 
be  used  as  a  proportionality  factor  for  the  entropy  generation  rate 
in  the  calculation  of  the  power  loss,  Eq.  (8),  has  remained  an  open 
question.  Lampinen  and  Wiksten  [43]  presented  a  solution  to  the 
classical  problem  of  how  to  determine  the  correct  temperature. 
They  explained  that  according  to  Holmberg  et  al.  [46],  if  the 
system  only  exchanges  heat  at  the  environmental  temperature, 
the  Gouy-Stodola  theorem  gives  the  correct  loss  of  power.  Most 
industrial  processes,  however,  exchange  heat  at  higher  or  lower 
temperatures  than  the  environmental  temperature  [46].  If  the 
environmental  temperature  is  then  used  to  calculate  the  loss  of 
power  in  these  cases,  the  Gouy-Stodola  theorem  does  not  give  the 
correct  loss  [46],  According  to  Holmberg  et  al.  [46],  using  the 
theory  of  Lampinen  and  Wiksten  [43] 


Ed  —  W  net.K—  W  ne[  —  TeffSgen 


where 


Teff  = 


h2~h2R 

S2—S2R 


(9) 

(10) 


Or,  if  the  specific  heat  capacity  and  pressure  remain  constant 


T2—T2r 

eff~\n(T2/T2R) 


(ID 


where,  for  a  system,  h2  is  the  specific  enthalpy  of  the  flow  at  the 
outlet  and  s2  is  the  specific  entropy  at  the  outlet  and  where  h2R 
and  s2R  represent  states  of  the  enthalpy  and  entropy  after  a 
reversible  process,  respectively. 

According  to  Bejan  [12],  it  is  important  to  stress  that  T0  in 
Eq.  (8)  is  the  temperature  of  those  areas  of  the  environment  that 
are  adequately  close  to  the  system  but  which  are  not  affected  by 
the  discharge  of  the  system.  According  to  Bejan  [12],  the  purpose 
of  such  a  modelling  decision  is  to  place  all  the  internal  and 
external  irreversibilities  related  to  the  physical  installation  inside 
the  system.  The  system  then  consists  of  the  installation  and  the 
surrounding  areas  which  are  affected  by  the  installation.  In  an 
actual  system,  rejected  heat  would  leave  the  installation  and  enter 
the  bordering  environmental  fluid  at  a  temperature  higher  than  T0. 
The  same  heat  would  also  cross  the  T0  boundary  further  down  the 
line  [12], 

Thus,  the  authors  of  this  paper  emphasise  that  the  choice  of  the 
control  volume  boundary  is  very  important  in  the  thermodynamic 
analysis  of  the  solar  thermal  Brayton  cycle.  The  authors  recom¬ 
mend  that  for  a  system,  Eq.  (8)  be  used.  Note  that  the  heat  losses 
from  components  should  then  be  calculated  as  the  heat  loss  up  to 
the  temperature,  T0.  When  calculating  lost  available  work  for  a 
single  component,  Eqs.  (10)  and  (11)  can  be  used  to  calculate  the 
boundary  temperature,  Tj  as  shown  in  Eqs.  (5)  and  (6). 


3.2.  The  sun  as  an  exergy  source 


Eqs.  (5)  and  (6).  The  following  literature  gives  a  review  thereof, 
by  highlighting  conflicts  and  new  perspectives. 

According  to  Izquierdo  Millan  et  al.  [47],  before  solar  radiation 
can  be  converted  into  mechanical  power,  it  is  first  converted  into 
heat  when  it  is  intercepted  by  an  absorbent  surface  located  on  the 
Earth.  This  occurs  with  a  great  loss  of  exergy.  The  exact  exergy  of 
solar  radiation  could  be  determined  with  spectral  measurement 
and  calculation  according  to  Petela  [48],  The  concept  of  solar 
exergy  maps  has  also  been  developed  byjoshi  et  al.  [49], 

According  to  Onyegegbu  and  Morhenne  [50],  the  expressions 
for  radiation  temperature  and  exergy  of  direct  beam  solar  radia¬ 
tion  are  well  documented.  The  expression  for  exergy  of  solar 
radiation  and  the  history  of  its  development  have  also  been 
investigated  and  documented  more  recently  by  various  authors 
[42,51-56],  Throughout  the  literature  [57-60]  the  absorbed 
solar  radiation  exergy  rate,  considering  the  Petela  theorem  [61], 
is  given  by 

^•HGtKGT) 

where  Ts  is  the  apparent  blackbody  temperature  of  the  sun  and  Q_ 
is  the  total  solar  power  input  crossing  a  control  volume  boundary. 
According  to  Eskin  [62],  the  expression  for  the  availability  of  solar 
radiation  (W/m2)  with  beam  and  diffuse  components  as  given  by 
Onyegegbu  and  Morhenne  [50]  is 

‘-'-HOM  '-m) 

where 

^  =  0.9562  +  0.2777  ln(l//c)  +  0.051  \k  (14) 

*  S 

The  expression  for  Ts/T*,  relating  the  effective  temperature  of 
diffuse  radiation  to  the  solar  temperature,  is  obtained  from 
Landsberg  and  Tonge  [63],  with  k,  the  dilution  factor  of  diffuse 
radiation,  less  than  0.1. 

According  to  Bejan  [8],  when  doing  an  exergy  analysis  on  a 
solar  thermal  system,  the  sun  can  be  considered  as  an  exergy-rich 
source  and  as  a  high-temperature  fuel.  Bejan  [44]  showed  that  the 
exergy  rate  of  the  sun  is  given  as 

£  =  Q*(l-p)  05) 

where  T*  is  the  apparent  sun’s  temperature  as  an  exergy  source 
and  Q  is  the  solar  heat  rate  crossing  a  control  volume  boundary. 
This  method  has  been  adopted  for  the  exergy  inflow  of  solar 
collectors  by  various  authors  [13,  33-36,  42,  64,  65],  Eq.  (15)  has 
also  been  used  to  describe  the  exergy  inflow  for  heliostat  surfaces 
[66,67],  The  abovementioned  authors  mostly  adopted  the  value  T*, 
as  suggested  by  Petela  [61],  being  approximately  equal  to  0.75TS 
[68,69]  (see  Eq.  (12)).  Ts  is  the  apparent  blackbody  temperature  of 
the  sun,  which  is  about  6000  K,  or  5762  K  [47,  70],  Therefore,  T*  is 
considered  to  be  close  to  4500  K  [68,  69], 

Farahat  et  al.  [58]  claimed  that  according  to  Najian  [71  ],  Eq.  (12) 
violates  the  second  law  of  thermodynamics  and  that  the  corrected 
equation  according  to  Torres-Reyes  et  al.  [72]  and  Najian  [71], 
assuming  the  sun  as  an  infinite  thermal  source,  is  Eq.  (16),  with  T* 
as  the  apparent  sun  temperature  and  equal  to  75%  of  blackbody 
temperature  of  the  sun  (4350  I<)  [68].  Also  note  the  inclusion  of  an 
optical  efficiency,  ij 


The  solar  input  for  the  solar  thermal  Brayton  cycle  can  be 
modelled  in  different  ways  in  an  objective  function  such  as 


E  =  ,,I 


(16) 
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According  to  Onyegegbu  and  Morhenne  [50],  the  expression  for  the 
exergy  flux  which  has  the  widest  acceptability  is  the  expression 

4-4’H©) 

The  selection  of  the  control  volume  around  the  analysed  system  and 
the  selection  of  Q  ,  or  the  term  in  front  of  the  bracket  in  Eqs.  (12) 
and  (15)-(17),  are  very  important.  Q  depends  on  which  boundaiy 
is  referred  to  in  the  analysis.  For  example,  Eqs.  (12)  and  (15)-(17) 
can  be  used  to  describe  the  incoming  exergy  from  the  sun  or  the 
incoming  exergy  from  a  reflector,  see  Fig.  2.  Q  can  thus  be  the 
beam  irradiance,  in  the  case  where  it  is  coming  from  a  reflector  for 
example  (hence  the  inclusion  of  the  optical  efficiency  in  Eq.  (16)),  or 
Q  can  be  the  global  solar  radiation  with  a  boundaiy  somewhere  in 
the  atmosphere  or  even  in  space.  For  example,  Xu  et  al.  [67] 
described  Q  as  the  incident  solar  radiation  so  that  Q  is  defined 
as  IA  where  I  is  the  amount  of  solar  radiation  received  per  unit  area, 
which  is  the  direct  normal  irradiation  (DN1).  The  available  solar 
radiation  per  unit  area  from  the  sun  varies  throughout  days,  months 
and  environmental  conditions.  The  best  database  for  solar  irradi¬ 
ance  would  be  the  long-term  measured  data  at  the  site  of  the 
proposed  solar  system  [73],  A  meteorological  database  can  also  be 
used  [74], 

According  to  [47,  70,  75,  76],  the  global  irradiance  or  the 
extraterrestrial  solar  radiation  normal  to  the  unit  area  is  the  solar 
constant  and  is  described  by 

S=foT4  (W/m2)  (18) 

where /is  the  sunlight  dilution  factor,  equal  to  2.16  x  10-5  on  the 
earth,  a  is  the  Stefan-Boltzmann  constant  (5.67  x  10~8  W/m2  K4), 
and  Ts  is  the  temperature  of  the  sun  (K).  According  to  De  Vos  [77], 
cited  by  [42,47,78],  solar  exergy  can  be  defined  as 

E=faT4-0~f)aTp4  (19) 

where  (1  —f)oTp  is  the  albedo  of  the  earth,  and  Tp  is  the  planet 
temperature.  The  term,  <tTC04,  can  also  be  subtracted  from  Eq.  (19) 
to  include  the  loss  of  exergy  due  to  the  temperature  of  the 
collector. 

According  to  [42,  70,  75],  Winter  et  al.  [79]  regard  the  exergy 
released  by  solar  irradiance  as 

£=/<rTs4(l-|^(l-0.281n/))  (20) 


Table  1 

Comparison  of  exergy  rate  calculated  with  different 
methods. 


Equation 

Solar  exergy  rate  (W/m2) 

(12) 

930.76 

(15) 

930.76 

(16) 

849.31 

(19)' 

904.84 

(20) 

985.54 

3.3.  Modelling  of  entropy  generation  mechanisms  in  the  solar 
thermal  Brayton  cycle 

In  a  solar  thermal  power  cycle  or  solar  thermal  convertor,  there 
will  always  be  entropy  generation  mechanisms.  These  mechan¬ 
isms  are  mostly  due  to  fluid  flow  and  heat  transfer.  Bejan  [8] 
provided  a  number  of  features  most  guilty  of  entropy  generation  of 
which  those  featuring  in  the  solar  thermal  Brayton  cycle  are:  heat 
transfer  across  a  non-zero  temperature  difference,  flow  with 
friction,  mixing,  filling  and  discharge,  compression  and  expansion. 
Bejan  [38]  also  presented  the  entropy  generation  mechanisms 
involved  with  the  conversion  of  solar  radiation  into  mechanical 
power  as  adiabatic  free  expansion,  the  transformation  of  mono¬ 
chromatic  radiation  into  blackbody  radiation,  scattering  and  net 
radiative  heat  transfer. 

Bejan  et  al.  [1]  considered  the  simultaneous  effect  of  heat 
transfer  and  fluid  friction  on  entropy  generation  and  showed  that 
they  have  a  tendency  to  contest  with  one  another  when  a 
thermodynamic  optimum  is  desired.  Entropy  generation  mechan¬ 
isms  have  also  been  described  by  Tolman  and  Fine  [45],  In  this 
section,  for  the  different  components  in  the  solar  thermal  Brayton 
cycle,  the  entropy  generation  mechanisms  are  considered  in  more 
detail  and  recommendations  are  made  on  how  these  entropy 
generation  mechanisms  can  be  put  into  equation  form  in  an 
objective  function  (see  Eqs.  (5)-(7)). 

3.3.1.  Solar  receiver 

According  to  Bejan  [8],  there  exist  three  main  features  that 
cause  irreversibilities  in  the  operation  of  any  solar  receiver:  heat 
exchange  between  the  sun  and  the  receiver,  heat  loss  from  the 
receiver  to  the  environment  and  the  internal  irreversibility  in  the 
receiver.  The  entropy  generation  rate  from  a  solar  receiver  can  be 
written  as  Eq.  (21)  [13,  38,  80],  where  the  solar  receiver  receives 
solar  radiation  at  the  rate,  Q 


Table  1  compares  the  methods  described  in  the  above  literature 
to  calculate  the  solar  exergy  rate.  It  is  assumed  that  Q  is  equal  to 
/=  1000  W/m2,  >/=0.85,  rt=Ts=5777  K  and  TP=300K.  Note  that 
Eqs.  (12)  and  (15)  give  the  same  result  for  all  practical  reasons. 

It  is  recommended  by  the  authors  that  Eq.  (16)  be  used  in  the 
analysis  of  the  solar  thermal  Brayton  cycle  with  T*  as  the  apparent 
sun  temperature  and  equal  to  75%  of  blackbody  temperature  of  the 
sun.  Also,  I  is  the  measured  DNI  and  >/  is  the  reflectivity  of  the 
reflector.  The  dish  reflector  is  thus  not  included  in  the  control 
volume  for  analysis. 


Fig.  2.  Control  volume  boundary  including  or  excluding  the  solar  reflector. 


Qo  ,  Qr_Q* 

To  +  Tr  T* 


(21) 


where  T*  is  the  apparent  sun  temperature  as  an  exergy  source,  T0  is 
the  ambient  temperature,  and  Tr  is  the  surface  temperature  of  the 
absorber  or  receiver  temperature.  For  a  non-isothermal  receiver 
and  neglecting  pressure  loss  [8,13,81] 


5  _Qo  ^  |  Tout  Q 

jgen  —  -f — h  niLp  III - — g 

1 0  l  in  I 


(22) 


The  entropy  generation  rate  in  the  receiver  of  an  open  and  direct 
solar  thermal  Brayton  cycle  is  given  in  Eq.  (23)  for  an  ideal  gas 
working  fluid,  where  the  effect  of  pressure  loss  is  included.  This 
equation  is  given  by  Le  Roux  et  al.  [33-35],  Note  that  Qo/To  in 
Eqs.  (21)-(23)  can  be  replaced  with  Q_ioss, receiver /T  b  where  Tb  is  the 
temperature  of  the  boundary  across  which  the  heat  is  being 
transferred.  With  this  replacement,  these  equations  can  be  used 
in  Eqs.  (5)— (7) 


q! 

T* 
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From  Eq.  (23),  it  is  shown  that  to  calculate  the  entropy  generation 
rate  of  a  solar  receiver,  the  net  absorbed  heat  rate,  heat  loss  rate, 
and  inlet  and  outlet  temperatures  and  pressures  are  required. 
These  values  will  depend  on  the  time  of  day,  month  and  year  as 
well  as  the  receiver  type  and  receiver  design.  Computer  software 
and  algorithms  are  available  to  compute  the  solar  heat  rate  on  a 
receiver  as  reflected  from  a  reflector  (see  for  example  [82]). 

Various  solar  receivers  can  be  used  for  the  solar  thermal 
Brayton  cycle.  These  include  tubular  cavity  receivers,  open  or 
closed  volumetric  receivers  and  particle  receivers.  Bertocchi  et  al. 

[83]  described  the  heating  of  air  temperatures  to  far  more  than 
1000  C,  using  a  high-temperature  solar  particle  receiver.  Kim  et  al. 

[84]  and  Miller  and  Koenigsdorff  [85]  also  described  the  features 
and  modelling  of  a  particle  receiver.  Kim  et  al.  [84]  found  many 
experimental  results  including  an  increase  in  the  opacity  of  the 
particle  curtain  with  mass  flow  rate  up  to  a  constant  value  near  the 
terminal  velocity.  Miller  and  Koenigsdorff  [85]  found  that  the  main 
loss  from  the  receiver  is  reflection  from  the  window,  followed  by 
emission  and  reflection  from  inside  the  receiver.  Heller  et  al.  [86] 
demonstrated  that  pressurised  (closed)  volumetric  receivers  are 
able  to  produce  air  of  1000  °C.  Pressurised  volumetric  receivers 
and  its  modelling  are  described  in  the  literature  [87-89],  Hischier 
et  al.  [87]  recommended  that  a  minimisation  of  the  cavity  wall 
thickness  in  relation  with  its  strength  to  withstand  the  operating 
pressures  should  be  performed.  Open  volumetric  receivers  were 
described  by  Avila-Marfn  [90]  and  Kami  et  al.  [91  ].  Kami  et  al.  [91  ] 
presented  a  “Porcupine”  volumetric  solar  receiver  and  its  ability  to 
withstand  a  solar  flux  of  up  to  about  4  MW/m2,  while  producing 
gas  exits  temperatures  of  up  to  940  °C.  This  receiver  consists  of  an 
array  of  pin-fins  or  elongated  heat  transfer  elements,  implanted  in 
a  base  plate. 

The  investigation  and  modelling  of  tubular  cavity  receivers  are 
available  from  the  literature  [92-94],  The  performance  of  different 
cavity  receivers  was  investigated  by  various  authors  [92,95,96], 
Shuai  et  al.  [92]  investigated  different  classical  cavity  geometries 
and  found  that  the  shape  of  the  cavity  (geometry)  has  a  significant 
effect  on  the  overall  distribution  of  the  radiation  flux  in  the  cavity 
receiver.  According  to  Shuai  et  al.  [92],  an  upside-down  pear  cavity 
might  be  an  appropriate  shape.  Prakash  et  al.  [95]  investigated 
heat  losses  from  a  cavity  receiver  at  different  inclination  angles, 
with  frontal  and  side  winds.  Reddy  and  Sendhil  Kumar  [96] 
numerically  compared  various  kinds  of  cavity  receivers  and  found 
that  their  modified  cavity  receiver  experienced  lower  convection 
heat  losses  than  other  receivers  and  suggested  that  their  modified 
cavity  receiver  may  be  favoured  in  a  solar  dish  collector  system. 
For  the  modified  cavity  receiver,  numerical  investigations  regard¬ 
ing  natural  convection  heat  losses  [97,98]  and  radiation  losses  [99] 
were  presented  for  the  modified  cavity  receiver.  Le  Roux  et  al.  [33— 
35]  described  the  modelling  of  a  modified  solar  cavity  receiver  in 
detail. 

When  a  dish  with  a  specific  diameter,  focal  length  and  rim 
angle  is  used  to  focus  sun  rays  onto  a  receiver,  the  net  rate  of  heat 
absorbed  by  the  working  fluid  in  the  receiver  depends  on  the 
aperture  diameter  of  the  receiver.  Due  to  the  sun's  rays  not  being 
truly  parallel  and  due  to  concentrator  errors,  the  reflected  rays 
from  the  dish  form  an  image  of  finite  size  centred  on  its  focal 
point.  The  aperture  area  of  the  solar  receiver  will  determine  the 
amount  of  intercepted  heat,  but  also  the  amount  of  heat  lost  due  to 
convection  and  radiation.  The  larger  the  cavity  aperture,  the  more 
heat  can  be  lost  due  to  convection  and  radiation  but  also,  the  more 
heat  can  be  intercepted  [4],  For  a  fixed  dish  concentrator  area,  the 
amount  of  heat  available  for  the  working  fluid,  which  is  the 
intercepted  heat  minus  the  heat  lost  due  to  radiation,  convection 
and  conduction,  is  a  function  of  the  cavity  aperture  diameter. 

Steinfeld  and  Schubnell  [94]  presented  a  semi-empirical 
method  to  determine  the  optimum  aperture  size  and  operating 


temperature  of  a  solar  cavity  receiver  for  maximum  energy 
conversion  efficiency.  They  found  that  for  inaccurate  concentra¬ 
tors,  with  Gaussian  flux  density  distribution  at  the  focal  plane,  the 
optimum  operating  temperature  varies  in  the  range  of 
800-1300  K.  Stine  and  Harrigan  [4]  presented  a  receiver  sizing 
algorithm  which  can  be  used  with  a  fixed  dish  concentrator  area  to 
get  the  intercepted  heat  for  a  specific  cavity  aperture  area. 

Using  the  abovementioned  literature,  the  authors  recommend 
that  Eq.  (23)  be  geometrised  (written  in  terms  of  receiver  geome¬ 
try)  for  a  specific  solar  receiver  design  with  specific  constraints,  so 
that  geometry  variables  can  be  used  in  an  optimisation  procedure. 


3.3.2.  Recuperator  or  heat  exchanger 

Heat  exchangers  play  an  important  role  in  the  solar  thermal 
Brayton  cycle  as  recuperator,  intercooler  and  in  indirect  systems 
where  the  working  fluid  is  not  cycled  directly  through  the  hot  side 
(receiver).  The  geometiy  variables  chosen  for  optimisation  in  a 
heat  exchanger  will  depend  on  the  design  and  type  of  heat 
exchanger,  since  many  different  types  of  heat  exchangers  exist.  A 
recuperator  is  used  in  the  Brayton  cycle  to  extract  the  heat  from 
the  turbine  outlet  and  transfer  it  to  the  cold  stream  before  it  is 
heated  by  the  heat  source. 

Entropy  generation  minimisation  has  been  utilised  in  various 
heat  exchanger  applications  for  optimisation  of  geometry 
[18,100,101].  These  optimisations  were  done  for  components 
individually  and  not  as  part  of  a  whole  system.  Lerou  et  al.  [101] 
showed  that  the  width,  height  and  length  of  flow  channels  can  be 
optimised  by  minimising  the  entropy  generation. 

According  to  Bejan  [8]  and  Bejan  et  al.  [1],  the  irreversibility  in 
a  heat  exchanger  is  the  sum  of  the  irreversibilities  of  each  of  the 
two  surfaces  of  the  heat  exchanger.  Two  factors,  temperature 
difference  and  frictional  pressure  drop,  are  responsible  for  irrever¬ 
sibilities  in  heat  exchangers  [8,102,103],  Yilmaz  et  al.  [103]  implied 
that  the  greatest  source  of  irreversibility  in  a  heat  exchanger 
comes  from  fluid  friction  in  the  form  of  pressure  drop. 

Exergy  analysis  has  been  applied  on  various  heat  exchanger 
types.  Sarangi  and  Chowdhury  [104]  expressed  the  entropy  gen¬ 
eration  in  a  counterflow  heat  exchanger  and  a  nearly  ideal  heat 
exchanger.  The  contribution  of  fluid  friction  to  the  entropy  gen¬ 
eration  was  neglected.  Adiabatic  ends  with  no  heat  loss  to  the 
surroundings  were  also  assumed  for  the  heat  exchanger. 

Exergetic  optimisation  has  been  done  for  tubular  heat  exchan¬ 
gers.  Cornelissen  and  Hirs  [105]  did  an  exergetic  optimisation  of  a 
balanced  water-to-water  counterflow  heat  exchanger  by  consider¬ 
ing  the  irreversibilities  due  to  pressure  drop,  due  to  temperature 
difference  between  the  hot  and  the  cold  stream  and  also  due  to  the 
production  and  construction  of  the  heat  exchanger.  Heat  loss  to 
the  environment  and  heat  resistance  of  the  tube  walls  were 
neglected.  Zimparov  [106]  included  the  effect  of  fluid  temperature 
variation  along  the  length  of  a  tubular  heat  exchanger. 

Second  law  analysis  and  optimisation  were  done  for  heat 
exchangers  with  ideal  gas  flow,  as  expected  in  a  solar  thermal 
Brayton  cycle  recuperator.  Hesselgreaves  [107]  considered  heat 
exchangers  with  zero  and  finite  pressure  drop.  When  looking  at 
zero  pressure  drop,  Hesselgreaves  [107]  included  balanced  coun¬ 
terflow,  flow  imbalance,  unbalanced  counterflow,  parallel  flow, 
condensing  on  one  side  and  evaporation  on  one  side.  Hessel¬ 
greaves  [107]  found  that,  for  zero  pressure  drop,  flow  imbalance 
increases  entropy  generation  and  that  it  is  advantageous  to  let  the 
highest  capacity  rate  stream  be  the  hot  stream.  Thermodynamic 
analyses  for  the  balanced  cross-flow  recuperative  plate-type  heat 
exchanger  with  unmixed  fluids  were  done  by  Ogulata  et  al.  [102], 
Ordonez  and  Bejan  [108]  did  a  numerical  optimisation  for  the 
parallel-plate  heat  exchanger. 
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The  entropy  generation  rate  in  the  recuperator  is  shown  in  Eq. 
(24)  for  an  ideal  gas  in  the  flow  channels  of  the  heat  exchanger. 
This  equation  or  a  similar  version  is  also  used  in  the  literature  by 
various  authors  [1,8,11,102,103,107,108]  and  is  recommended  by 
the  authors  of  this  paper 

Sgen  =  rriCp o  ln(T2/Ti)  +  rhCpo  ln(T4/T3)+rhcp0  In  [(P2/Pi)(1“k,/,t] 

+  mcp0  In  [(P4/P3)(1-k)/fc]  (24) 


The  fluid  going  from  position  1  to  position  2  is  the  cold  stream  and 
the  fluid  going  from  position  3  to  position  4  is  the  hot  stream  in 
the  heat  exchanger  as  shown  in  Fig.  3. 

According  to  Ordonez  and  Bejan  [108],  entropy  is  also  gener¬ 
ated  due  to  the  discharge  at  the  recuperator  in  an  open  cycle— 
hence  Eq.  (25).  The  equation  to  be  used  in  an  analysis  depends  on 
the  definition  of  the  boundaries  of  the  control  volumes  of  the 
recuperator  and  the  system 


Sgen  =  rhCpo  ln(T2/Ti)  +  mCp0  \n(T4/T3)+mCp0 


T2-To 

To 


+  mcp o  In  [(P2/Ptf~k)/k]  +mcp0  In  [(P4/P3)(1^k)/fc] 


+ 


Q/o: 


T„ 


(25) 


Eq.  (25)  should  not  be  used  when  the  entropy  generation  due  to 
the  exhaust  of  the  open  system  (warm  air  outlet)  is  already  added 
as  a  separate  entropy  generation  mechanism.  The  authors  recom¬ 
mend  that  Eq.  (24)  is  used  in  conjunction  with  an  objective 
function  such  as  Eqs.  (5)-(7). 

Different  designs  are  available  for  the  Brayton  cycle  recupera¬ 
tor.  In  solar  applications,  the  recuperator  is  often  designed  as 
integral  to  the  micro-turbine.  To  model  the  recuperator  entropy 
generation  (Eq.  (24)),  the  inlet  and  outlet  temperatures  and 
pressures,  and  the  heat  loss  from  the  recuperator  are  required. 
The  heat  loss  to  the  environment  from  the  surface  of  the 
recuperator  can  be  a  significant  factor  and  it  is  recommended  that 
it  is  included  in  solar  thermal  Brayton  cycle  analysis.  Hesselgreaves 
[107],  Ogulata  et  al.  [102]  and  Ordonez  and  Bejan  [108]  suggested 
that  the  e—NTU  (effectiveness -number  of  transfer  units)  method, 
based  on  the  second  law  of  thermodynamics,  can  be  used  to  get 
the  outlet  temperatures  and  the  total  heat  transfer  rate  from  the 
hot  fluid  to  the  cold  fluid.  This  is  shown  in  the  literature  by 
Haugwitz  [109]  and  Le  Roux  et  al.  [33-35], 

In  many  counterflow  heat  exchangers,  the  heat  flow  through 
the  material  in  the  longitudinal  direction  is  neglected  in  determin¬ 
ing  the  temperature  profile  over  the  heat  exchanger.  Eq.  (26) 
([101])  is  a  dimensionless  parameter  that  can  be  used  to  see  if 
longitudinal  conduction  can  be  neglected  or  not.  k  is  the  thermal 
conductivity  of  the  heat  exchanger  material  and  A  is  the  cross- 
sectional  area.  The  longitudinal  conduction  cannot  be  neglected  if 
Abh  >  10~2  [101] 


2bh  = 


\<A 


Lrhcn 


(26) 


A  few  examples  of  recuperator  designs  and  the  modelling  thereof 
are  given:  heat  exchangers  are  required  to  be  efficient,  safe, 
economical,  simple  and  convenient  [103],  Heat  transfer  and  pres¬ 
sure  losses  as  well  as  the  optimisation  of  price,  weight  and  size 
should  be  considered  while  designing  the  heat  exchanger  [102], 
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Fig.  3.  Heat  exchanger  with  cold  stream  (1-2)  and  hot  stream  (3-4). 


According  to  Bejan  [8],  heat  exchanger  irreversibilities  can  be 
decreased  by  slowing  down  the  fluid  which  is  travelling  through 
the  heat  exchanger.  Also,  for  a  fixed-area  heat  exchanger,  the 
irreversibility  can  be  reduced  by  allocating  the  area  correctly  [8], 

Kreith  and  Kreider  [110]  and  Hesselgreaves  [107]  suggested 
that  counterflow  heat  exchangers  should  be  used  in  solar  thermal 
power  systems  and  that  parallel-flow  heat  exchangers  should  be 
avoided.  According  to  Bejan  [8],  counterflow  heat  exchangers  are 
often  used  in  recuperative  heating  associated  with  the  Brayton 
cycle.  It  was  also  suggested  by  Shah  [111]  that  counterflow  plate- 
type  heat  exchangers  can  be  used  as  compact  recuperators  with 
micro-turbines.  Shah  [111]  gave  design  criteria  for  micro-turbine 
recuperators.  Criteria  such  as  high  performance  with  minimum 
cost,  high  exchanger  effectiveness,  compactness,  40,000  h  opera¬ 
tion  life  without  maintenance  and  low  pressure  loss  (  <  5%)  were 
given.  According  to  Shah  [111],  these  criteria  translate  into  a  thin 
foil  primary  surface  recuperator  where  flow  passages  are  formed 
with  stamping,  folding  and  welding  side  edges  by  an  automated 
operation.  Such  a  recuperator  was  discussed  by  McDonald  [  112,113], 

Pra  et  al.  [114]  described  printed  circuit  technology  and  plate 
fin  technology  for  recuperators  while  Tsai  and  Wang  [115]  inves¬ 
tigated  the  design  and  analysis  of  a  Swiss-Roll  recuperator.  Le 
Roux  et  al.  [33-35]  described  the  modelling  of  a  plate-type 
recuperator.  Traverso  and  Massardo  [116]  discussed  the  furnace- 
brazed  plate-fin  type  and  the  welded  primary  surface  type 
recuperators.  Ultiainen  and  Sunden  [117]  reviewed  a  recuperator 
with  cross  corrugated  or  chevron  pattern  heat  transfer  surface. 
Different  recuperator  designs  [118]  and  the  effects  of  recuperator 
channel  geometries  [119]  are  available  from  the  literature.  Yilmaz 
et  al.  [103]  presented  the  literature  on  optimisation  based  on 
entropy  generation  numbers,  for  the  following  heat  exchangers: 
balanced  and  unbalanced  counterflow  heat  exchangers,  cross-flow 
heat  exchangers,  external  flow  heat  exchange  (like  fins),  two- 
phase  flow  heat  exchangers,  regenerative  heat  exchangers,  plate- 
type  heat  exchangers  and  shell-and-tube  heat  exchangers.  Heat 
exchangers  with  a  liquid  can  be  modelled  with  the  help  of 
references  [1,  8,  103,  106,  108,  120-122],  Such  heat  exchangers 
can  be  found  in  applications  such  as  exhaust  gas  heat  recovery  or 
intercooling. 

Various  attempts  have  been  made  at  entropy  generation  mini¬ 
misation  for  a  specific  component  or  a  specific  entropy  genera¬ 
tion  mechanism,  by  defining  the  entropy  generation  equation 
as  a  function  of  geometry  variables  [1,8,120-122].  Eq.  (24)  for 
the  recuperator  or  heat  exchanger  entropy  generation  rate  in  a 
solar  thermal  Brayton  cycle  can  be  geometrised  similarly  using  the 
literature  above. 

Using  the  abovementioned  literature,  the  authors  recommend 
that  Eq.  (24)  be  geometrised  (written  in  terms  of  recuperator 
geometry)  for  a  specific  heat  exchanger  design,  so  that  geometry 
variables  can  be  used  in  an  objective  function  such  as  Eqs.  (5)— (7). 


3.3.3.  Compressor  and  turbine 

Eqs.  (27)  and  (28)  show  the  equations  for  the  entropy  genera¬ 
tion  rates  for  the  compressor  and  turbine  in  the  case  of  an  ideal 
gas.  These  equations  were  also  given  byjubeh  [11]  and  Bejan  [8], 
Eqs.  (27)  and  (28)  are  formulated  in  terms  of  the  pressures  and 
temperatures,  which  could  be  described  by  the  isentropic  efficien¬ 
cies.  The  definitions  of  these  efficiencies  are  commonly  available 
from  Weston  [123],  Dixon  [124]  and  Haugwitz  [109]  for  example. 
The  geometries  of  compressors  and  turbines  are  often  complex 
and  the  authors  of  this  paper  do  not  recommend  that  the  entropy 
generation  rate  of  the  compressor  and  turbine  be  written  in  terms 
of  geometry  variables  of  the  turbine  or  compressor.  The  calculation 
of  entropy  generation  rate  for  these  components  at  different 
temperatures  and  pressures  is  however  important  in  the 
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derivation  of  the  objective  function.  The  pressures  and  tempera¬ 
tures  can  rather  be  written  in  terms  of  the  geometry  variables  of 
other  components 


By  assuming  that  the  net  power  output  (generator  output  power) 
would  be  equal  to  the  acceleration  power  at  no  load,  it  can  be 
shown  that 


mcp0i_e  In (Te/Ti)-mR  ln(Pe/P,)+^£ 

(27) 

mcpoi-e  ln(Te/Tj)-mR  In  (Pe/P)i+^|^ 

(28) 

In  this  paper,  the  modelling  of  radial  flow  turbines  is  discussed  in 
more  detail  since  it  is  mostly  used  in  small-scale  power  cycles. 


3.33.1.  Compressor.  The  compressor  entropy  generation  rate  can 
be  calculated  with  Eq.  (27).  The  compressor  pressure  ratio,  mass 
flow  rate  and  isentropic  efficiencies  are  shown  in  a  compressor 
map,  given  by  the  compressor  manufacturer.  These  maps  are 
mostly  fitted  with  functions  to  make  the  modelling  easier  as 
shown  by  Haugwitz  [109]  and  Frei  [125],  or  fitted  using  software 
such  as  GTPower  or  TCMAP  or  as  shown  by  Westin  [126].  Methods 
to  model  the  compressor  based  on  its  geometry  are  also  available 
[127],  Wahlstrom  and  Eriksson  [128]  also  presented  methods  of 
compressor  modelling. 

Haugwitz  [109]  presented  a  method  by  Gustafsson  [129]  wherein 
continuous  ellipsoid  curves  were  fitted  to  speed  curves  on  the 
compressor  map.  According  to  Haugwitz  [109],  the  parameters  in 
the  ellipsoid  equation 


can  be  adjusted  to  fit  any  speed  curve  on  a  typical  compressor  map. 
For  parameterisation,  the  compressor  efficiency  map  can  be  mod¬ 
elled  as  parabolic  degradation  curves  [129],  The  amount  of  degrada¬ 
tion  or  curvature  (d)  is  fitted  based  on  data  to  get  parameterisation 
for  all  speeds.  According  to  Gustafsson  [129]  the  following  equation  is 
valid  for  one  speed  [109]: 

7c  =  dc,  max- — nrmax  e^-)  (30) 

where  ;/c  max  is  the  maximum  efficiency  at  that  speed  and  rhmax  eff is 
the  corresponding  mass  flow  for  this  maximum  efficiency.  Other 
compressor  modelling  options  are  described  in  the  literature  by 
various  authors  [126,130-133], 


3.3.32.  Turbine.  Turbine  maps  are  not  as  available  as  compressor 
maps  and  even  if  they  are  available,  the  turbine  efficiency  is  not 
always  shown  on  the  map,  but  rather,  a  maximum  efficiency  is 
given.  This  often  makes  turbine  modelling  difficult.  According  to 
Watson  and  Janota  [134]  and  experimental  results  by  Westin 
[126],  Haugwitz  [109]  and  Zhuge  et  al.  [127],  it  can  be  assumed 
that  the  turbine  efficiency  is  a  function  of  the  pressure  ratio  and 
the  speed  of  the  shaft  and  that  a  maximum  efficiency  exists  for 
each  speed  line.  According  to  Reuter  et  al.  [135],  the  instantaneous 
isentropic  turbine  efficiency  can  be  calculated  with 


7t(f)  =  bmWr/t.s  = 


W  c(  t)  +  W  acc(  t ) 

wts(t) 


(31) 


where  Wt,s(t)  is  the  instantaneous  isentropic  turbine  power  and 
the  acceleration  power  can  be  determined  with  Reuter  et  al.  [135] 
and  Westin  [126]  as 


Wacc(t)  = 


2 

JrotorNrotor(t) 


dNrotor(t) 

dt 


(32) 


Westin  [126]  suggested  that  the  polar  moment  of  inertia  can  be 
determined  with  the  trifilar  method  as 


(33) 


rnCpT7iit(1—rt(k~1Vk)—^-mCpTj(rc<k~1)/k—1)  =  (J£)  )rotor NrcoAt)^^ 

(34) 


This  equation  can  be  used  to  determine  the  turbine  efficiency  if  J  is 
available  and  if  a  time  dependant  modelling  is  used.  Another 
method  to  determine  turbine  efficiency,  using  the  blade  speed 
ratio,  is  discussed  in  the  following  paragraph. 

According  to  Westin  [126],  Renberg  [136],  Watson  and  Janota 
[134]  and  Moraal  and  Kolmanovsky  [131],  steady  flow  theory 
states  that  if  the  data  in  an  original  turbine  map,  showing  a  curve 
on  axes  of  pressure  ratio  vs.  mass  flow  rate,  is  plotted  as  efficiency 
vs.  blade  speed  ratio  (H/Cs  or  BSR )  all  data  points  would  fall  onto  a 
parabolic  curve,  with  a  definite  peak.  This  peak  should  be  posi¬ 
tioned  at  U/Cs= 0.707  if  the  ratio  of  expansion  between  turbine 
rotor  and  housing  (reaction)  is  50%  which  seems  to  be  common  for 
vaneless  radial  turbines.  This  is  also  found  in  the  doctoral  thesis  of 
Shaaban  [137  who  suggested  that  the  blade  speed  ratio  is  the 
most  commonly  used  parameter  to  define  the  turbine  efficiency. 
This  phenomenon  seen  when  plotting  efficiency  vs.  BSR  is  also 
found  and  shown  by  turbocharger  testing  and  research  [138],  It  is 
debatable  how  the  turbine  works  under  unsteady  flow— Westin 
[126]  showed  that  the  blade  speed  ratio  at  maximum  efficiency 
could  stay  away  from  the  standard  value  of  0.707,  at  non-steady 
conditions.  According  to  Watson  and  Janota  [134]  the  problem 
can  be  treated  as  quasi-steady  flow,  meaning  that  the  turbine 
performs  under  non-steady  flow  in  the  same  way  as  it  would  if 
those  sudden  flow  conditions  were  steady.  Evidence  on  the  error 
introduced  by  this  assumption  is  not  entirely  consistent  and  it  is 
commonly  believed  that  the  error  will  not  exceed  5%  [134], 
Capobianco  and  Marelli  [139]  showed  that  the  errors  are  due  to 
this  assumption.  The  blade  speed  ratio  can  be  determined  as 
[126,128,132] 


BSR  = 


V(D\ 

rlz) 


[2hi(1-rtV-kVk)]'/2 


(35) 


so  that  when  it  is  assumed  that  the  maximum  efficiency  is  found 
at  BSR=  0.707,  rt  can  be  found  for  a  specific  speed  (N).  The  points  of 
maximum  turbine  efficiency  for  each  speed  line  can  be  plotted  on 
the  given  turbine  map  (mass  flow  rate  vs.  pressure  ratio).  This 
maximum  efficiency  is  also  given  on  the  turbine  map.  When 
considering  experimental  results  by  [109,126,127]  it  is  assumed 
that  the  turbine  efficiency  can  be  presented,  with  pressure  ratio  on 
the  x-axis,  by  plotting  curves  for  each  speed  line  which  goes 
through  the  point  of  maximum  efficiency.  The  efficiency  can  be 
modelled  as  a  parabolic  function  of  the  blade  speed  ratio  [125, 127, 
128,  140] 


7f  —  It,  max 


BSR— 0.6\  2\ 

.  0.6  J  J 


(36) 


According  to  Guzzella  and  Onder  [133],  simplified  closed-form 
descriptions  of  the  turbine's  efficiency  can  be  approximated  by 


7r (BSR)  —  }jr  max 


/  2  BSR 
\BSR0pt 


(37) 


where,  in  automotive  applications,  typical  values  are  rit.max  ~  0.65- 
0.75  and  BSRopt  x  0.55-0.65.  Thus,  Eqs.  (36)  and  (37)  are  essen¬ 
tially  the  same. 

From  the  above  it  is  concluded  that  from  a  basic  turbine  map 
the  turbine  efficiency  can  be  calculated  as  a  function  of  pressure 
ratio  and  speed  with  the  use  of  the  blade  speed  ratio.  The  authors 
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of  this  paper  propose  that  this  method  will  aid  the  analyst  to  a 
great  extent  in  determining  the  entropy  generation  rate  of  the 
turbine. 


3.3.4.  Radiator 

The  radiator  is  used  at  the  cold  side  of  a  closed  solar  thermal 
Brayton  cycle.  The  entropy  generation  rate  depends  on  the  type  of 
heat  exchanger  used  for  the  radiator.  Eq.  (38)  is  given  for  the 
radiator  if  fins  are  used  to  cool  down  the  working  fluid  (being  an 
ideal  gas)  with  a  natural  external  flow  (being  air) 

Sgen  =  y^+ifiCpo  ln(j^)  ln(pr)  (38) 


where  the  subscript,  B,  refers  to  the  base  where  the  fins  are 
attached.  According  to  Bejan  [8],  Eq.  (38)  can  also  be  written  as 


Q.B.sngiaB-Too)  ,  FdV . 


+- 


+  JhCpo  ln(  7re)-mR  In  fj 


(39) 


where  it  is  assumed  that  Qi0ss, radiator  =  Qb  =  NQ.B,sngi  and  N  is  the 
number  of  fins. 

There  is  competitiveness  between  heat  transfer  and  fluid 
mechanics  terms  to  get  an  optimal  radiator  size  with  minimum 
entropy  generation  [1,8,17], 


3.3.5.  Other  entropy  generation  mechanisms 

Other  entropy  generation  mechanisms  occurring  in  the  solar 
thermal  Brayton  cycle  include  entropy  generation  due  to  mixing 
(refer  to  [8]),  filters,  pipes,  pipe  bends  and  storage.  The  entropy 
generation  due  to  these  mechanisms  when  assuming  an  ideal  gas 
flow  can  be  calculated  with 

■Sgen  =  ^Y^+rhcpo  In  (^\  —mR  In  (jy^J  (40) 

The  modelling  of  fluid  flow  in  straight  pipes,  bent  pipes  and  flow 
splits  is  discussed  in  the  literature  [109,136].  Other  typical  flow 
restrictions  can  be  modelled  as  shown  by  Frei  [125],  Phase  change 
material  has  been  modelled  by  Wu  et  al.  [141]. 

Often,  the  bearings  of  turbo-machines  are  cooled  using  oil  or 
water.  The  entropy  generated  from  such  a  system  can  be  described 
as  shown  in  Eq.  (41)  and  included  in  an  objective  function.  The 
cooling  flow  loop  can  also  be  modelled  separately  or  as  the  heat 
loss  in  Eqs.  (27)  and  (28) 

Sgen  =  %^  (41) 

1  b 

When  a  solid  or  a  liquid  crosses  the  control  volume  boundary  of 
the  system,  the  entropy  change  can  be  modelled  with  Eq.  (42) 
instead  of  Eq.  (4) 

se-Si  =  c  In  ^  (42) 


4.  Conclusions  and  recommendations 


thermodynamics,  including  geometry  into  optimisation  and  the 
importance  of  the  selection  of  the  control  volume  boundary,  was 
emphasised.  The  modelling  of  various  components  to  construct  an 
objective  function  for  optimisation  was  discussed.  Such  an  equa¬ 
tion  can  be  applied  for  example  to  optimise  the  component 
geometries  in  a  steady-state  solar  thermal  Brayton  cycle  such  that 
the  system  produces  maximum  net  power  output. 

Various  studies  have  stressed  the  significance  of  the  optimisa¬ 
tion  of  the  global  performance  of  a  system,  by  minimising  the  sum 
of  the  irreversibilities  from  all  the  different  components  or  pro¬ 
cesses  of  the  system  (distributing  the  entropy  generation  rate 
through  the  system  by  optimally  sizing  the  hardware,  instead  of 
optimising  components  independently)  [1,17,18,38,108,121,122], 
For  future  work  it  is  recommended  that  geometry  optimisation 
studies  be  done  on  especially  intercoolers  and  reheaters  (a  second 
solar  receiver)  in  the  solar  thermal  Brayton  cycle,  as  the  inclusion  of 
these  components  can  largely  increase  the  efficiency  of  such  a 
system.  The  authors  also  recommend  that  constraints  on  the 
objective  function,  such  as  component  size,  cost  and  weight,  should 
be  incorporated  depending  on  the  component  types  chosen  in  an 
analysis. 

In  this  paper  it  was  shown  that  from  a  basic  turbine  map, 
turbine  efficiency  can  be  calculated  as  a  function  of  pressure  ratio 
and  speed  with  the  use  of  the  blade  speed  ratio.  It  is  proposed  that 
this  method  will  aid  the  analyst  to  a  great  extent  in  modelling  and 
determining  the  entropy  generation  rate  of  the  turbine  in  the  solar 
thermal  Brayton  cycle.  The  inclusion  of  turbine  and  compressor 
geometry  variables  for  optimisation  in  the  objective  function  can 
be  considered  for  future  work.  The  authors  referred  to  different 
types  of  receivers  and  heat  exchangers  in  this  paper  and  recom¬ 
mend  that  these  components  can  be  modelled  and  optimised  in 
terms  of  total  system  entropy  generation  minimisation  for  com¬ 
parison  purposes.  Entropy  generation  in  the  transient  state  was 
not  shown  in  this  study  and  can  be  recommended  for  future  work. 
Other  entropy  generation  mechanisms  in  the  solar  thermal  Bray¬ 
ton  cycle,  not  mentioned  in  this  paper,  and  the  modelling  thereof 
can  also  be  recommended  for  future  work. 

The  work  in  this  paper  is  applicable  to  solar  thermal  studies  in 
general  but  focused  mainly  on  the  small-scale  solar  thermal 
Brayton  cycle.  With  this  paper,  the  authors  hope  to  contribute  in 
making  solar  thermal  Brayton  systems  successful  and  applicable  in 
future  power  generation. 
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In  this  paper,  new  perspectives,  conflicts  and  recommendations 
were  highlighted  regarding  the  derivation  of  an  objective  function 
which  includes  various  geometries  from  different  components  of  a 
solar  thermal  Brayton  cycle.  Optimisation  studies  of  the  Brayton 
cycle  as  relevant  to  the  solar  thermal  Brayton  cycle  were  reviewed. 
The  influencing  factors  to  optimisation  were  identified.  The 
method  of  total  entropy  generation  minimisation  was  highlighted 
as  it  allows  heat  transfer  and  fluid  flow  terms  to  be  available  for 
optimisation  in  a  single  equation,  so  that  geometry  variables  of 
different  components  of  a  solar  thermal  Brayton  cycle  can  be 
optimised  simultaneously.  The  importance  of  the  second  law  of 
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